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ABSTRACT 


The perception hns been that the equatorial Jets observed n JupUer and 
Saturn cannot be understood In the framework of a eonally symmei rlc circula- 
tion; that some kind of eddy process equivalent to upgradlent d f fusion wlte 
negative viscosity Is required to acceli'rate the motions. We a rive at a 
different conclusion when vm examine this problem with spectral models 
describing a viscous i luld. Near the equator where the Corloll^ force 
vanishes, the momentivn Ktlance Is established between horizontal and vertical 
diffusion, which, a ptlocl, does not Impose constraints on the «'lrectlon or 
magnitude of the zonal winds. Solai ra llatlon absorbed at low Mtltudes Is a 
major force In driving large scale muclons with air rising near the equator 
and falling at high latitudes. In the upper leg of the uerldloi il cell, 
angular nKuaentum Is n llstrlbuted so that the atmosphere tends 'o subrotate 
(or corotate) at low latitudes and superrotate at high latitudes. In the 
lower leg, however, the process Is reversed and produces a tend* ncy for the 
equatorial region to » -iperrotate. The outcome depends on the energy budget 
which Is closely coupl *d to the momentum budget through the ther.aal wind 
equation; a pressure (temperature) maximum Is required to sustain equatorial 
superrotation. Such u condition arises In regions wlilch are convect I vely 
unstable and the temperature lapse rate Is superadlabatlc. It should arise 
In the tropospheres oi Jupiter and Saturn; planetary energy from the Interior 
Is carried to higher altitudes where radiation to space becomes Important. 
Upward equatorial motions In the direct and Indirect circulations (Ferrel- 
Thomson type) Imposed by Insolation can then "trap" dynamic enerity for 
equatorial heating to sustain the superrotation of the equ.itorlal region. Duo 
to vertical momentm diffusion and convection the dynamic signatures from the 
troposphere and lower stratosphere are "mixed" and produce a complex 
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Interference pattern In the zonal circulation near the tropopaua< . We 
attribute the large dUferencea between the clrculatlona on Jupiter and Saturn 
to differences In alr.e, gravitational accele rat Ion, temperature, eddy 
dlffualvlty, and tropospheric stability* 
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INTRODUCTION 


AtmoHpherlc super rotation is observed on virtually every pl.met In our 
solar system (e.g* Rln);-Hele, 1964; Conrath et al., 1973; Smith rt al., 

1979a, b; Alexander, l'i(i2; Schubert et al. , 1980). In each case, the planet's 
rotation axis is nearlv perptendicular to the orbital plane around the Sun, the 
eq'jatorial region Is hi'ated preferentially and its pressure decriasts toward 
higher latitudes. Applying thj concepts of geostrophlc and cyclost rophlc 

equilibria, reasonably sm.ill latitudinal pressure gradients (ord>T of 

1 0 

i if' o-» can account tor \ wide range of superrotatioii rati'S that are 
observed. On the basis of this evidence, we are led to believe 'hat thenniil 
forcing (e.g. Lorenz, 1967; Leovy, 1973; Glerasch, 19753 Is a major 11 not the 
prime cjindldate for driving atmospheric superrotation. 

In an earlier pap< r (Nayr and Harris, 1981 ) we lud consldei'ed the 
angular momentum budget of a thermally driven zonally symmetric circulation 
and provided an elementary description of the processes leading to and 
maintaining the rigid shell component of supcrrotatlon. It was thown that 
depending on the transport properties, an atmosphere can go through two 
different modes. For 'arge viscosity It behaves almost like a sulld body, 
frlctlonally coupled t the planet. The angular momentum fur superrotation Is 
then "derived" from th decrease In the moment of inertia due to mass 
redistribution by the Meridional circulation; during the "spin up time" the 
atmosphere transfers angular momentum to the planet. For lower <nd more 
realistic values of viscosity there Is an Important additional degree of 
freedom In the vertical momentum distribution. Mass can more freely return 
from the poles In the lower leg of the meridional cell and the c'lange In the 
moment of Inertia Is not as large. This return flow together with the 
Increased tendency toward geostrophy causes subrotation In the lower regions 


of the atmosphere during the Initial phase of the spin up process. Through 
viscous shear the planet applies a torque to the atmosphere which propagates 
upwards and in time supplies the angular momentum for large superrotation. 

Thus in steady state, uiider thermal forcing, a rigid shell component of 
superrotation eventually prevails at all altitudes. 

A longstanding and remaining problem is differential rotation. Out- 
standing examples are the banded wind fields and the very large well defined 
equatorial jets on Jupiter (e.g. Smith et al., 1979a, b) and Saturn (Smith et 
al., 1981) which are not understood. 

The general perception is that equatorial sUiK rrotatlon cannot develop in 
a thermally driven zoniilly symmetric circulation te.g.. Hide, 1969, 1971; 
Gierasch, 1975; Held an! Hou, 1980). Large scale eddy processes ire therefore 
invoked to explain the ibservatlons. For Jupiter it is estimated (Ingersoll, 
19ol) that eddies ranging In size from 100 to 1000 km contain enough energy to 
accelerate the prevailing zonal circulation within a short period of one 
month, giving credence to the notion that , somehow, "order comes out of 
chaos." 

We suggest in this piaper, contrary to current theory, that the equatorial 
Je';s on Jupltar and Sat irn can be understood In the framework of a thermally 
driven zonally symmetric circulation. First, an empirical model Is presented 
for Jupiter's zonal win! field observed by the Voyager spacecraft. The power 
spectrum in terms of vector spherical harmonics conveys a very simple and 
beautiful order, In which momentum can be’ perceived to cascade from lower to 
higher order modes. Second, the thec-em of equatorial subrotation (e.g., 

Hide, 1969; Held and Hou, 1980) Is crltlrally examined and an alternative 
interpretation is presented. Near the c<|uator where the Coriolis force 
vanishes, the momentum balance is established between horizontal and vertical 
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diffusion, which, a pilorl, does not impose constraints on the direction or 
magnitude of the sonal winds. Thus, the problem is reduced to explain how 
momentum and energy are transported Into the equatorial region? Third, a 
spectral model and the concept of mode coupling are Introduced ti describe the 
direct and Indirect circulations driven by the absorption of sol. r radiation 
at low latitudes. Fourth, the momentum and energy budgets are qralitatlvely 
discussed for the meridional circulation driven by solar differential 
heating. Relative to the rigid shell component of superrotation which 
receives most of Its angular momentum from the planet during spin up, angular 
momentum Is redistributed In the upper leg of the merldionil cell so that the 
atmosphere tends to suhrotate or corotate at low latitudes and superrotate at 
high latitudes. In the lower leg, however, the process Is reversed and a 
tendency for the equatorial region to superrotate Is produced, '■/hether this 
actually occurs depends crucially on the energy budget which Is closely 
coupled to the motnentiim budget through the thermal wind equation, a pressure 
(temperature) maximum Is required to sustain equatorial superrot itlon. The 
necessary condition for this should arise In the tropospheres of Jupiter and 
Saturn, where small scale convection must transport energy from the Interior 
upwards, and therefore the temperature lapse rate should be superadlabatlc. 
Upward equatorial motions In the direct and Indirect circulations (of the 
Ferre 1-Thomson type), Induced by Insolation, can then "trap" or "focus" 
dynamic energy for equiitorlal heating which Is required to sustain equatorial 
superrotation. Fifth, using paramet erlzatlons In terras of Newtonian cooling 
and Raleigh friction, simple analytical solutions for the spectral model are 
presented which elucidate some basic properties of the zonal circulation In 
different regimes of the atmosphere and under different planetary conditions. 
Sixth, allowing for molecular and eddy transport processes. Including vertical 
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and horlEontal nomentum diffusion in the meridional and zonal foi ce balanceSi 
a numerical solution is presented vhich describes equatorial sup« rrotation in 
a zonally symmetr^v' circulation driven by solar differential heating. Angular 
momentum is thereby allowed to cascade from lower order modes wh< re most of 
the energy is absorbed to higher order modes. The integration carried out 
from the lower troposphere (25 bar pressure) up into the exosphei e. 


OBSERVATIONS 

Spectral analyses of the zonal wind field on Jupiter observed by the 
Voyager Spacecraft (G. Hunt, private communication) are shown in Figure 1. 

They are based on empirical representations of the data In terms of vector 
spherical harmonica which also serve as basis functions for our . heoretlcal 
model (Mayr and Harris, 1981). The amplitude-^ In the lowe. part of Figure 1 
are the equatorial val es of the solenoldal harmonics (Morse and 

Feschbach, 1953). Positive and negative values contribute to su 'errotatlon 
and subrotation at the equator respectively. In the upper part »f Figure 1 
analytical representat ons of the Voyager I data from both hemls|>heres are 
presented for synthese up to order L-12 and 34. Included are only the 
harmonics which are symmetrical with respect to the equator (even wave 
number i ). 

The fundamental h rmi>nlc (H*2) or che rigid shell component of super- 
rotation dominates. Higher order modes, i>4, are at least a faci or of two 
smaller arJ are apparei'tly in a separate class of their own. Af' er the abrupt 
decrease between i-2 and l»4 the amplitudes decrease gradually and monoton- 
Ically toward higher wave numbers. Near i»18 the sign changes, ind the 
spectrum up to t«30 contributes to equatorial subrotation. We note that the 
spectra for the Voyager I and II observations are nearly Identic 1. 
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Anticipating the theoretical r«!8ult9 discussed later we >ffer the 
follo%rlng interpretation. The rigid shell component of superrot itirn (1-2) Is 
primarily driven by differencial heating from solar radiation wb ;rc most of 
the energy is deposited in the lowest order harmonics (e.g., Vol land and Hayr, 
1972a). We refer to ir as "direct circulation". That component receives most 
of its angular momentum from the planet or the planetary interior during spin 
up (Mayr and Harris, 1^81 ). In contrast, the solar input into the higher 
order modes is relatively small and "direct forcing" is probably not 
important. This part )f the circulation spectrum is excited by aode coupling 
in which momentum (and energy), originating primarily in the solir driven 
rigid shell component >f superrotation, cascade from lower to higher order 
modes. We refer to it is "indirect circulation". 

For vector spherii'al harmonics, the Coriolis force induces node coupling; 
the range of momentum • ransfer In t space being proportional to the planetary 
rotation rate. By comparison with the global scale circulation on Venus, a 
slowly rotating planet, the zonal velocity fields on Jupiter and Saturn are 
therefore confined into n^irrow latitude bands. From the synthesis in Figure 1 
one can see that the 1 )wer order harmonics (L-12) primarily contribute to tie 
equatorial jet, suppor’ Ing our interpretation that solar radiation absorbed In 
the lowest order harm< ilcs is driving these motions. We attribute equatorial 
superrotation to momentum advection In the lower leg of the merl Ilona 1 cell 
and energy convection into the region oelow the cloud top, where, presumably, 
the temperature lapse rate is supendiabatlc. The negative (equatrrlal) 
velocities in the power spectrum around <-22 are largely responsible for the 
positive Jet near 25*^ latitude. We believe that this feature is the signature 
of poleward momentum and energy transport in the upper leg of the meridional 
cell. 1716 effect is not (or to a much lesser extent) seen on Saturn (Smith et 
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(il>, 1931), Indicating that the clouiU form at a lower altitude c >ep Inside 
the (unstable) convection region, consistent with the lower temp'^ratures 
observed on this planet (Hanel et al. 1981). 


DIFFUSIOt. ; LARGE SCALE CIRC UL ATION 

Assuming |U| << , the conservation of angular momentum with respect 

to the planetary rotation axis can be expressed as 


^•(V M ) 
o 




AM, 


+ (Fj) 


( 1 ) 


where 


M 

o 

AM 


■ (U, V, W), wind vector with zonal, meridional and \ertlcal 

components 

2 ’ 

•’ (i»^r p^cos'0, "planetary" angular momentum of the atmosphere 

- Urp cosO , atmospheric angular momentum 
ngP© 32 jj 

■ -X ( S'), momentum diffusion (viscous flow) due tr 

m gg2 rcose 

meridional shear 


p^, global averat'.e of obss density, 
H, density scale height, 
m, HBSS 


n , n , vlscosltv coefficients In the vertical or horizontal directions, 
r 0 

(i)^, planetary au;alar velocity 

r, radial dlstan' e from the center of tl>e planet 

0 latitude 

From the continuity of mass flow, ^ hand side of (I) 

reduces to (VM ) ■ rVp slri0cos0 which vanishes at the equator. 
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Following Held aid Hou (1980), the aasuatptton Is mad-* that Fq 0> Since 

3 U 

U-0 the planetary aurface and yp ■ 0 at the top rf the atnos )here, it 
follows that 

U (r) - o (equitor) (2) 

at all altitudes. 

A parcel of air -novlng In a Hadley cell from the equ-itor t >woi ; the pole 
is accelerated under conservation of angular momentum giving rl <e to super- 
rotation. Nowhere along the trajectory can the total angular mxnentum, 

Mj. - + AM, exceed the planetary angular momentum at the equitor (Eq. 2). 

This is an upper limit considering viscous dissipation, and it follows 


U < 


b) rsln‘0 

_P 

cosO 


(3) 


Equation (3) is due to Hide (1969) and expresses a fundamental theorem of 
atmosphere dynamics. Consiacrlng the constraint (Eqs. 2 and 3) on the zonally 
symmetric circulation. Hide (1969), Held and Hou (1980), Tngersill et al. 
(1979) and others con.-luded that eddy processes, equivalent to -ipgradlent 
diffusion with negatl .'o horizontal viscosity, must provide the ingular 
momentum for the equatorial accelerations on Jupiter and Saturn. 

We question here the general validity of Hide's theorem (E]. 3). By 
neglecting hori«:untal momentum diffusion (F^ ■> 0) in (1), only the Coriolis 
force remains to balance vertical diffusion, and that balance breaks down near 
the equator, producing the trivial result that the atmosphere nist corotate 
(2). But given the physical reality of horizontal diffusion, all that can he 
said about the equator Is that horizontal Influx or out flux of momentum must 
be balanced by vertical outflux or Influx of momentum, respect 1 v'e ly. A 
priori, under this force balance, there is no constraint on the direction or 


9 


.-niignltude of Che zonal wind velocity; the aCeoephere la p«>niltt 'd either to 
suhrotaCe or auperrotate at Che eqitaCort A good example la the obaerved wind 
field on Jupiter (Flg-ire 2). At latltudea below about 7®, angular mooentun la 
dlffualng downgradlent toward the equator where It can be carrl *d through 
vertical shears to lower and or higher altltiides. Thla being t'ue, the 
queatloii reduces to h >w angular monentm la transported Into or out of the 
’’equatorial region" (Figure 2)? 

We assume chat solar radiation preferentially absorbed at low latltudea 
supplies the energy f >r the thermally driven meridional circulation with air 
rising near the equator and falling at higher latitudes* Under conservation 
of mass, more planetary angular rooroentum Is carried upward by V'rtlcal motions 
at low latitudes than Is descending at high latltudea. Ttils gl >1 m 1 Inbalance 
Is conpensated by downward diffusion or viscous flow (e.g., Leo/y, 1973; 
(llcra-’ch, 1975) which requires the atmosphere to assume a global average cr 
rigid shell component of superrotat Ion all the way down to the surface (Figure 
3). The angular mons ntum la supplied primarily from the ,'lanet during spin 
up. Given that basic component, the question Is what are the c uitacCerlstlcs 
of differential rotaUor related to differential (equatorial) heating? 

MOMENTUM AND EKER~Y B VLANCg 

In answering the -e qtiestlons we present first an idealized and heuristic 
Interpretation of equatorial superrotat 1 on. At altitudes above the center of 
the meridional clrculitlon, where the heat Input per mass Is relatlvlely 
large, angular momentim Is effectively advected poxeward, excep' in tlie 
"equatorial region" where the meridional wlnJa and the Corlolla force gf» to 
zero and the horlrontil momentum t rinapurt muat be carried by h»rlzontal 
diffusion. Thla causes the atmospheric rotation rate, relative to the rigid 
shell component of superrotation, to decrease (minus) at low latitudes and 
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Increiise (plu.') at high lAtlCudes (Figure 4). The equatorial r glon therefore 
tends to corotate or subrotate, consistent with Hide's th< oren Eq. 3). Relow 
the center of the nerldional circulation the process is reverse. . Angular 
noacr*'un is advected from high to low latitudes so that, relati'e to the rigid 
shell conponent of sui'errotation, the rotation rate tends to de> rease at high 
latitudes and increase at low latitudes (Figure 4). Again, in the vicinity of 
the equator, it is understood that the momentum transport Is ta! en over by 
horizontal diffusion. 

However, this picture is not unique as demonstrated in the upper part of 
Figure 3. Three different velocity distributions are 3hown whlih can balance, 
through vertical diffusion, the momontum advoction into and out of ti>e 
equatorial region in the upper and lower legs of the meridional cell. At the 
planetary surface or deep Inside the atmosphere the velocity is assumed to he 
zero. For case (b) the zonal velocities are positive (superrot. lion) and 
negative (subrotation) at lower and higher altitudes, respectively, 
representing the scentrlo described In Flgur.* 4. But lor (o) and (c) the 
atmosphere superrotates and subrotates at all altitudes. 

Whethei the redl^ t rlbution of angular momentira by the meridional circula- 
tion actually leada to superrotat icn or subrotation crucially depends on ^he 
energy budget which lx closely coupled to the momentm budget through the 
thermal wind equation. Applying perturbation theory, we consider the conser- 
vatl.on of energy In simplified forr 

.2 3T 

(a - X 1^) AT + c w ( ° > D - AQ (4) 

r ,2 p 3i 


where 



specific heat per luss at constant pressure 
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a, coefficient of radiative cooling, 
thermal conductivity per unit mass 
W, vertical velocity 

Tg, average temperature 

At, temperature perturbation 

r, adiabatic lapse rate 

AQ, energy Input per unit mass 

For tlie Indirect circulation (AQ ■ 0), the temperature perturbation is 
determined by energy convection due to vertical motions 


AT 


3T 

-c w(^ + r) 
p 9r 


3r 


(5) 


Under "local" condlti >ns It is reasonable to assume >c — - < 0. 

3T^ 3r* 

lorraally, with si atlc stability ^ 0» the upward uotlors (W > 0) 

near the equator caus.- the temperature to decrease, AT < 0. Tl'rough the 
thernvil wind equation, the resulting pressure decrease tends to produce 
subrotation, and cast (c) in Figure 5 describes the moment nm budget. This 
condition applies to he atmosphere at higher altitudes. 

However, when the temperature lapse rate Is superadinbatlc , 


3T 

(yy- + r ) < 0, the Indirect meridional circulation Is capable of "trapping" 
dynamic energy conveci.ed upwards fror: lower altitudes (Figure ?). TTie 
temperature then Incr *ase8, AT > 0, and the accompanying pressure maximum 
permits superrotation In the equatorial region. Such a condition arises In 
regions of the atmosphere which are cliaracterl?ed by global scale energy 
diffusion cn higher altitudes wliere radiation to space Is Important. On Venus 
this Is observed below the visible cloud top, between 50 and 60 kra. It should 
also arise In the tropospheres of Jupiter and Saturn where energy from the 
planetary Interior Is carried upwards by small scale conv''Mor To a much 
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lesser extent , it nsy also arise In the troposphere on Ea^th. 

In summary, we emphasize some subtle points which ari* Important to under- 
stand equatorial superrutatlon: 

(a) Balancing vertical momentum diffusion with Coriolis force, when 
applied to the eqi.atur, led to an unphyslcal constraint on the tonal wind 
velocity (Eq. 2). Wl'h the reality of horizontal diffusion, this process i!uy 
or may not be Important In various regions of the atmosphere, lut It Is 
certain that horlzont il diffusion must he Important near the equator. The 
vanishing Coriolis force on Its own cannot provide a balance fo- vertical 
diffusion which la known to play a major part In the globril momentum budget. 
With the Inclusion of horlzoutal diffusion, there Is no constraint on the 
direction or magnltuJo of the equatorial zonal wind velocity. 

(b) Under the constraint (2), the tacit assumptions were mide that the 
meridional motions ca ry angular Twraontiun from the equator toward the poles 
(only), and that In t ils process angular momentum Is (always) 1 >st through 
diffusion; the consequence of the argument w:is Hide's theorem (1). At higher 
altitudes where the s>eclflc heat Input Is relatively large, it Is Indeed 
reasonable to assume i hat the meridional motions are directed nvay from the 
equatorial source and that dissipative processes due to dlffuslm would tend 
to dampen the zonal u(nd velocities. However, In the lower leg of the 
meridional clrculatloi where the return flow occurs, this process is 
reversed. The motion, are directed toward the equator. Angular momentum 
diffusing downward# fiom higher altitudes can be picked up and transported to 
lower latitudes which hMuld tend to accelerate the equatorial region. A 
factor contributing t>- this problem may have been the neglect o'' viscous 
moment urn transfer In the meridional force Ivtlance. Under tliat simplifying 
assumption the meridional velocltleb cannot be forced to vanish, as they 
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should, at the planetaty surface. Thus an ex'esslve amount of M.e momentum 
return flow la confined to the lowest altitudes where the ronal -'elocltles 
must vanish. 

(c) Finally, the I’lcture of a fluid belnt; spun up by polewa d motions 
under conservation of .lagular momentum was Incomplete for anothei reason, the 
most Important one per'iips. Vfhlle vertical diffusion was consld> red with the 
Inclusion of dissipation, momentum and energy convection by the ^•ertlcal 
motions In the morldloi'al circulation did not play a j>art In Hide's theorem 
(3). Yet we know that this transport process Is Important In supplying the 
atmosphere with angulai momentum from the vast rt^servul r ol the planet; 
computer simulations o< the spin up In which the average or rtgli shell 
component of siiperrotai Ion Is eslabllsh-'d cloirlv dem>>nstraLe ^ Importance 
^Mayr and Harris, 1*1(11 ). Within the atmosphiTe thif angular mo lent uo Is then 
redistributed by the 1 direct clrcul.it Ion, .m.l the upward notion, at the 
equator play again an mp irtant part. In regions ot the at nosph -re, below 
radiating cloud-or ha?. -1 • 'rs where the temp<>rature lapse rate an be 
supcradlabat l c, energy convection from below Increases the terape ature (and 
pressure). This Is tin u provldee tlie hisls, energetically, for iccelerating 
the equatorial region Irectly, or Indirectly due to vertical expansion and 
dlCfiu.lon wlilch can ca rv angular momentum to higher altitudes. 

MODEL AMD MODK COl^LIN 

To describe dlfftMentlal rot.itlon wi' adopt a spectral model (Mayr an 1 
Volland, 1972; Ilirrls nd Mayr, 197S; Mayr et a1., 1973; r«*ferre 1 t»> as M(l) . 
Its prop''rLles relevan tor our .inaWsts are bilelly summ.i 1 xet' • (1) 

Advectlon and convection of energv, nuss .nid nomiiit'irD associated with the 
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large scale zonally symmetric circulation are considered. (2) /ertlcal and 
horizontal diffusion >f energy and momentum (In both the ional tnd meridional 
force balances) are 1 icluded. (3) Perturbation theory Is applied relative to 
a globally uniform at losphere which corotatos with the underlyl ig planet; the 
physical variables su h as temperature, density and wind vector are expanded 
In terms of spherical harmonics (with the zonal wave number bel ig zero) which 
effectively separate latitude from altitude variations. (4) Using tridlagonal 
block elimination, th-- resulting set of dlfterence equations ar* solved from 
the planetary surface (or atmospheric Interior) up to the exosp lere. (5) At 
the lower boundary, the vertical, meridional and zonal velocltlis as well as 
the density and tempe-ature variations are set to zero. For th » upper 
boundary conditions, i he dominance of molecular diffusion requi *es that the 
vertical gradients of temperature and horizontal winds vanish. 

No separable elg>-n functions exist to describe the latltud; dependence 
for the prevailing atiospherlc circulation. Hough modes are appropriate eigen 
functions only wnen d sslpatlve processes due to viscosity can !>e neglected 
(Chapman and Llndzen, 1968), and spherical harmonics are appropriate only when 
the Coriolis force Is unimportant (Volland and Mayr, 1972b). Neither of these 
conditions is satisfied, and mode coupling must be considered. 

For a set of basis functions, chosen to be spherical harmonics, this 
Implies that an exter lai source In one spherical harmonic "mode” not only 
excites an atmospherli' perturbation In that particular mode; It also excites, 
to various degrees, all the other modes. Given a latitudinal distribution of 
the source, the structure of the atmospheric response changes depending on the 
importance of mode coupling, and one of the principle determining factors In 
that Is the planetary rotation rate (Coriolis force). 

Since Individual modes are coupled together, their definition should. In 
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principle, not effect the solution <is long as the set of basis functions Is 
complete and the expansion Includes a sufficiently large number of terns. For 
laatheioatlcal convenience we chose ”P^ modes" of order t each di fined by a set 
of (vector) spherical himonlcs (Pj^, Sj, Morse and Feschl »ch, 1953) 


AQ \ 

AT I 

Alogp I 
W / 


I (-1)^'^ X, (r) P. (0) 
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> V «/•> 

V - I (-1)*'^ (r) 

, A . . . 


- I (- 1 ) 

t-2,A... 


(i/2-H ) 


Vi Vl 


(7) 

( 8 ) 


wliore Alogp Is the l> ;arlthm of the pressure ratio (disturbed lo undisturbed). 
Foi even t values th« circulation Is symmetrical with respect to the 
equttor. The ampllti les In (6), (7) and (8^ are so defined tlv»t , arbitrarily 
close to the equator, positive amplitudes (r), (rl, (t)) i?lve 

positive cont rlbutlottn to the physical variables such 
as AT (r, 01, V (r, ( > and U (r, 01. 

With this ansnt • the equations of enetgv, mass and momentum conservation 
take the form 
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where g Is the gravlt tlonal acceleration anJ 


(4') 


2 /4'(4'+2)(4'+l)^ 
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For brevity, we ire omitting here a number of tenna Involv ng vertical 
and horizontal energy and oomentvm diffusion. In the lowest or ler mode, for 
example, the horlzont I momentum flux should vanish. The complete expressions 
are considered In our numerical an;)lysls. 

The equations (9 through (lA) describe the mode In the context of a 
larger system schemat cally Illustrated in Figure 6. On the right hand side 
of (12) and (13) are he momentum sources due to mode coupling which arise 
from the Coriolis for e and Involve the adjacent modes. 

To simplify the llscusslon, we assume that the atmospheric circulation Is 
forced only with the lowest order heat source 


- -C.2 (r) ?2 

AQj - 0, for 4 > 2, 


(15) 


representing differential heating by solar radiation which Is preferentially 
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absorbed at low latitudes. 


In an earlier paper (Mayr and Harris, 1981 ) we discussed s >me 
characteristics of the fundamental ?2 mode or the rigid shell coiponent of 
superrotation which Is directly driven by the solar Input (15). "rtils mode Is 
unique In that Its glo’>al angular momentum and moment of Inertia do not 
vanish. During spin up, that angular momentum la primarily supplied by the 
planetary torque at the surface and Is transported to higher altitudes through 
viscous shears and the meridional circulation. In contrast to the fundamental 
mode, the global angul ir momentum and moment of Inertia of the higher order 
modes P. (i>2), describing differential rotation, are Identically rero, and the 
problem Is no longer where the angular m«)menium for superrotation comes from 
but how It Is redistributed within the atmosphere? 

Liter on a numerl •al soluton ol tlie differential equations (9) through 
(13) will be presented In which we consider In self consistent form the mode 
(momentum) coupling or glnatlng from the lowest order heat source (15). 


QUALITATIVE ANALYSIS 

The vertical structure of the atmospheric circulation is pu ramet erl red 
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where the sign of Is determined by the sign of V|^ In the low>r leg of tne 
meridional clrculstlo i. h la a measure for the vertical extent of the 
circulation, and Is positive* 

Since thermal forcing by solar radiation Is only considers 1 In the lowest 
order mode (15), It 1> reasonable to assume, that. In the proce is of cascading 
momentum to higher or ier modes, the zonal and meridional wind onponents would 
tend to decrease 

"t-i - “t+i ■ 

h-2 - V, - V, (1 - 6) (21) 

The parameter 0 < 0 < 1 vanishes as mode coupling ceases to be Important. 

Substituting (Ih) through (21) and Ignoring momentum diffusion In (12\ 
the equations (9) thr >ugh (14) take the fotm 

3T 

<a-*c, 7 c +r) |V, I . q, ( 22 ) 

r 

^ h AT^ - mm(Ujj_j- ~ um(l-B)U^_j (23) 

o 

(t(^ -Ig + C^> Ut.,- »» (Vt _2 - V,) - «(1-B) V, tor t . 2 

- -tarn Vjj for t - 2 

The equations (22) and (24) describe the conservation of energy and angular 
momentum respectively, and the thermal wind equation (23) describes the 
coupling between energy and angular momentum. To first approximation (23) 
does not Involve diffusion processes and Is therefore more “rigid" than (22) 
or (24), where or can be either positive or negative. Normally, unlesa 
specified otherwise, It Is reasonable to assume that local conditions prevail 


19 


and 


/ } >0. (25) 

For the lowest order lode, t"2,the heat source q^>0 cause.4 the r.empera- 
ture to Increase at the equator AT2>0< Equation (23) then yields thu rigid 
shell component of superroi aMon, U^>0, which Is consistent with the momentum 
balance (24) for meridional winds flowing avay from the equator (V2<3) in the 
upper leg of the clrculatl m* In the lower leg of the circulation the dlrec* 
tlon of the meridional wlnls Is reversed. Put the rigid shell component Is 
unique In that, In steady itate, It cannot exchange angular momentum with the 
planet. To satisfy this constraint, the vertical flux of momentum cist 
approach zero near the low<*r boundary requiring that the sign of the curvature 
Cl7) must change (C^<0). :onslstent with the energy and momentum balances, a 
rigid shell component of i- iperrotatlon thus prevails all the way down to the 
planet (Figure 7a). 

For 1>2 and q^j^ - 0, t le situation Is different. In such an "Indirect 
circulation" the motions nre driven through momentum coupling between adjacent 
modes, and the temperatun- variations arise from dynamic heating. 

Moreover, for higher order modes there Is no constraint on the viscous stress 
at the lower boundary, ~ — ) * 0. We consider the indirect circulation to 

be forced by solar (equat rial) heating :'.n the lowest order mode, q^>0. At 
low latitudes the merldloivil winds In the lower leg ut the cell are then 
directed toward the equator, and the po'jltlve sign applies In Equation (22). 

The single most Important factor determining the dynamics Is the 
3T 

term + T In the energy balance which Is also a measure for the 

stability of the atmosphere. With S^>0 the atmosphere Is convert Ively stable. 
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the temperature lapse rate la subadlabattc and» in the global average » heat la 
diffusing downwards* Conversely, with S^<0 the atmosphere is unstable, the 

lapse rate Is superadlabatlr and thermal energy la diffusing upwards. 

3T 

(r) ■ + r > 0: This condition holds, when, above r, the 

atmosphere absorbs more radiative energy than It emits. It usually applies to 
the more tenuous regions at higher altitudes, far away from the heat 
dissipation by the planet, where the solar radiation Is less attenuat :d and 
colllslonal excitation Is li*ss effective. With q^-0, the upward motions at 
the equator then cause the temperature to decrease, AT^^<0 (22), and (23) 
requires that the zonal wind component also de:reases (U£_j< 0) relative to 
the rigid shell component of superrotation. This Is consistent with poleward 
transport of angular momentum In the upper leg of the circulation (24). The 
return flow of angular momeotun toward the equator In the lower leg o*^ the 
cell, however, tends to Increase the zonal velocity so that Its curvature must 
change sign (C^<0), fulfilling force balance (24). Such a sceairlo Is 
Illustrated In Figure 7b. 

S^(r) ■ + r < 0: This condition holds when, above r, the atmosphere 

emits more radiative energy than It absorbs. It usually applies to tlie denser 
regions at lower altitudes, c'^'se to the heat dissipation by the plan>‘t, where 
the solar radiation is more attenuated and colllslonal excitation Is -oore 
effective. With q^-0, the i.pward motions at the equator cause the tenperature 
then to Increaae, 4T^^> 0 (2.‘), and (23) requires that the zonal wind >'omponent 
also Increases 0) relative to the rigid shell component of 

superrotation. This Is con-lstent with the equatorward return flow of angular 
momentum In the lower leg of the circulation (24). In the upper leg of the 
cell, however, angular momentum Is transported atray from the equator which 
tends to decrease the zonal velocity to the exte.it that Its curvature must 
change sign (C^<0), fulfilling force balance (24). Such a scenario is 
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illustrated In Figure 7c. 

'3rt(r)<0 (lower altltudt^c) and S^(r)>0 (higher altitudes) ! Under 
realistic conditions it is reasonable to essuae that planetary .-itaoepheres arc 
convectively stable with su'>adlabatic lapse rates at higher altitudes and that 
they would tend to be conveotirely unstable with superadiabat^c lapse rates at 
lower altitudes. For such i scenario* the circulation is then \ combination 
of the two earlier exaraples and its structure near the equator is Illustrated 
in Figure 7d. At lower altitudes* below the tropopause (S^<0)* the upward 
notions supply energy (22) ^hich is carried toward the eq<iator in the lower 
leg of the meridional cell. As a result the temperature increases (&T^>0) and 
forces the zonal wind compo lent to Increase (U,_,>0) relative to the rigid 
shell component of euperrot ition. This is consistent with momentum transport 
In the lower leg of the mer dional circulation which is balanced (24) by a 
divergence in the vertical lomentin diffusion (C^>0). Above the center of the 
meridional cell and above tie tropopause (S^>0)* energy is carried away from 
the equator and the tempera’:ure decreases (ATj<0), forcing the zonal velocity 
to decrease relative to the rigid shell component of superrotation. 

The poleward transport of mimentum is again balanced by a convergence in the 
vertical momentm diffusion (C^>0). 

The picture portrayed here is somewhat simplified. In reality* the 
pressure increase due to thermal expansion in the troposphere extends upwards 
into the region which is convectively stable and thereby transports angular 
momentum to higher altitudes. Moreover, vertical diffusion of energy and 
momentum also contribute to couple the different regions. The result is that 
a complex Interference pattern can develop such as we see it in the fovlan 
ri rculation. 

For a net emitted flux of radiation, F^, and given a vertical eddy 
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dltiuHion ccMi'ricicnt ch«> "stability" at tiio t>rii>tauro Uvul It 
estliaaLrd 


3T 

So' <77^ + ■■> 


kT P 
o o 

c p X* 

P o 


(2f»> 


On Jupiter, - l.bxlO^ergi/ca^aec (Orton and Ingoraoll, 1976). Assuming K • 
10^ and Tq « 170^t ^ derive a value of ■ -10 deg/cu .‘i: the I b.^r 
pressure level which proves significant for the zonal circulation. Comp>ared 
to Jupiter, ^Mttuiu lb lutther away from the Sun and Che net catwiro tlu< of 
redldUlon In bis'ller (R. H^nel, private conmunlcatlcn). But, Its troposphere 
appears to be less turbuli t. Both effects would tend to compensate each 
other so that the superadlabatlc lapse rate, on Saturn may be comparable to 
(nsL cf Jupiter, un Venus, superadlabatlc lapse rates have been Inferred for 
Che region between ii) to 6C* km below the visible cloud top (Schubert et al., 
ivbu) which, we believe. Is Important to understand the equatorial super- 
cotation on that planet (M yr and Harris, 1981). 

L atitude S tructure ; To understand the nature of mode coupling iff acting 
the latitudinal structure of the circulation, we discuss here a simplified 
urtalytlcal wolucloo for tho equations (22) through (24). This will also 
elucidate some of the impo tant elements that distinguish the atmospheric 
circulations oi> Jupiter, S-iturn and Venus. 

for q^ • 0, we consld<;r the homogeneous equations (22) through (24) and 
oiicaiu the expr«bsion 






rujm(l-0) / H 

P O (_^y^ 4C ) 




(27) 
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Ignoring radiative cooling and aaaualng tg**- 

the coupling coefficient can be eatloatad 


B - 



( 28 ) 


For email t values, ne^tr the fundamental mode where solar energy is aep.. sited, 
mAi^ntum coupling la mi>st effective; the "slippage** bettreen the merlui.inal 
velocity components of adjacent modes Is small, and the value for u is cljac 
to one. Conversely, with 8*^0, mode coupling becomes email and a clu.'actet- 
iKCic uhve number R can be defined, which la a measure for tne LHiigt «ver 
wnlch momentum caacade»< 



We consider the lower leg of the meridional circulation wltli 

Vj^<0 and S^<0. Adopting h«l00 km, Tq« 170° and 8^—10“^ for Jupiter, 

360 ^ 

uotaiiis R-28 and a characteristic latitude structure of X ■ — ■ U**, in 

reasonable agreement with observations (Smith et al., 1979a, b). For lM.curn, 
the gravitational acceleration as well as the planetary radius and t.lacion 
rate are sniller than for Jupiter. If we assume that on both planets w* 
t aa;.cratures and the s>iperadlabatlc temperature lapse rates are about t..e same 
at the cloud level, this yields R - 14 or X - 24^, consistent with t .«a Voyager 
measurements (Smith et al., 1981). The atmosphere of Venus la m>ite couareased 
and we assume h ■ 10 kii. Adapting an atmospheric rotation (lerloc tu i h« order 
of 10 days and S^ • we obtain R • 4 or X - 90°; the clrculailuU «a 

global scale character, In substantial agreement with obusr vat lout t*>b><lieit et 
el., 1980). 

The Equations (22) through (24) may serve to estimate t'..*- 
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mod« or tha rigid iholl coapononc of tupocrotatlon 


V, - 
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( 30 ) 
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Conoldorlng sod* coupling In (24) duo Co V^_ 2 > taking V 2 tram (31), tho 
•oluClon for Cho hlghor ordor lodes (q^*C, i>2) can Chon bo obtained in tho 
fora of olmplo rocurolon rolaclanuhlpo glvon bolow. 
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(34) 



^■(l-B) Hr ,, 
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(35) 


Tho vortical velocltj lo given by (19), and the coupling coofflclent 0 lo 

defined vltb (2V) o^ (28) In olapllflod fora. 

* 0 * •“ 

A idJiP rigoioue apptoach to tho lu .o...ogci>oouii (q2^0) equatlono (22) 

through (24), trlLhout the need of a coupling coofflcl'sct 0, would require a 

• • • • . 

n'isprical ocluLlop thre alloiro for ootaontivi tranafor over o wide lange of 


2: 


modes. Tne order of this system would be limited by horizontal diffusion 
which damnens the amplitudes at rates proportional to — tc„ 


r 

The analytical model (27) through (35) proves useful to gain some Insight 

Into the properties of atmospheric circulations. We present here an example 

-4 

representative of the conditions on Jupiter. Adopting u ■ 1.8x10 , Sq ■ 

±10”^, m ■ 3.3x10“^^, H ■ 100 km and K » 10^ and Ignoring radiative cooling, 
results from the model are shown In Figure 8. Presented are the power spectra 
for the computed amplitudes of the zonal velocity, at the equator and 

their synthesis describing the latitudinal variations. It Is assumed that 
under negative and positive stabilities, S^, the meridional winds are directed 
toward the equator (Vjj^>0) and the poles (Vjj<0) at lower and higher altitudes 
respectively. For positive values of and C^, the expression under the 
radical (27) is then positive. 

Compared to the fundamental mode or the rlgla shell component of super- 
rotation (Ujj^>l) .which Is driven primarily by the solar Input (q^ >0), the 
Indirect circulations driven by the momentum coupling with the lowest order 
mode are distinctly different. At lower altitudes wliere Sq< 0, the Indirect 
circulation contributes to equatorial superrotation. The amplitudes 

starting vlth £ » 4, are significantly smaller than that of the 
fundamental mode and fall off gradually toward higher wave numbers (shorter 
wave lengths). This Is In substantial agreement with the observations 
for £ < (Flgart i). In tha eyntuesis cf this spectrum a narrow equatorial 
jet Is reproduced. 

At higher altitudes with the direct clrcilatlon Is again driven by 

differential solwr heating whlrn pr'^duces s jperrotatlon (U.>0). However, the 
Indirect ■ Irculatlon Is now characterized by energy and moiflentum transport 
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totmrd tha polas* Starting with tha tonal valoclty co ^‘'nanta ara naga- 
tlva and contrlbuta to aquatorlal subrotatlon. Synthaalslng thla power 
apactnia, tha moult la that tha ataoapham naarly corotatoa at tha equator 
and auparrotataa at aldlatltudea. Tha Inclusion of horlsontal diffusion and, 
partially, tha truncation In our apactcal analysis account for aooentua 
transfar toward the equator at low latitudes which sustains a snail 
suparrotatlon rate. 

Due to '-‘srtlcal convection and diffusion, the above discussed signatures 
from the upper and lower regions of the atnosphera are mixed, and we suggest 
that this Is responsible for the complex clrculatlou part'.'n observed on 
Jupiter. The observed negative equatorial velocity components around 1"22 In 
Figure 1 may originate from a spectrua of the kind shown with tha theoretical 
results In the lower part of Figure S. 

With the Large rotation rate and large planetary radius of Jupiter, the 
second term In the denominator of (30) dominates. Assuming O"0, the 
meridional velocity from (30) Is V 2 *q 2 > and the rigid shell component of 
super rotation from (31) is U^*q 2 /K* Considering (28), (33) and (34), a trend 
analysis for the Indirect circulation with t>2 gives 



where Y Is an expression, usually positive, which Involves atmospheric and 
planetary parameters. W:>.th Increasing a major component of the power 
spectrin Is shifted toward lower wave numbers. 
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NUMERICAL RESULTS 


We are demons era ting vlth this paper that In principle and contrary to 
current theories equatlonal superrotation can he understood In the framework 
of a tonally symmetric circulation driven by solar differential heating. ’Ms 
suggests that the process Is Important for the large scale atmospheric motions 
observed on Jupiter and Saturn. To determine the extent that this Is true 
will require a nonlinear analysis with general circulation models which 
account In self consistent form for the solar Input and convective as well as 
radiative energy transfer. Hie results will depend on the vertical and 
horizontal eddy dlffuslvltles and on the PranJtl number, and transport of 
latent energy la probably Important. 

Our discussion here Is restricted to convey a general understanding of 
the mechanism. No attempt Is made to simulate a particular planetary 
scenario. 

For Illustrative purposes a global average atmosphere la adopted which 
rotates at a rate of 1 .3xl0“^sec“^ In a gravltat loruil field with g-2400 and Is 
representative of the conditions on Jupiter. The Aita are based on the 
Infrared measurements from the Voyager spacecraft (lUnel et al., 1979a, b). 

They are shown In Figure 9. Below the 1 bar pressure level, the temperature 
Is extrapolated to lower altitudes assuming a (negative) stability • -S x 
10”^ whose absolute value decreases toward lower altitudes satisfying flux 
continuity. The lower boundary la qhosen to be at a pressure of 25 bar. 

Above the 2 mb level we assume the temperature to bo constant. The upper 
boundary Is chosen to be at a pressure of 5x10”^^ b.ir where dissipative 
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processes due to molecular diffusion dominate. 

For the energy source driving the circulation, Aq^Cr), we chose a height 
distribution shown in Figure 9 which was adopted from the work of Wallace et 
al. , 1974. The specific heat input increases with height and accounts for the 
temperature inversion above the tropopause. Radiative cooling is considered 
in the form of a Newtonian coefficient a-6.5 x 10 ^p^ (ergs deg"^ sec“^gra“^) 
provided by B. J. Conrath (private communication). Vertical and horizontal 
eddy diffusion coefficients of 2.5x10^ and 2.5x10^^, respectively, are chosen 
at the lower boundary and are assumed to increase with height at a rate 
Inversely proportional to the square root of the density (Chapman and Lindzen, 
1970). 

We have used two different schemes to obtain solutions of the 
differential equations (9) through (13): 

(1) Iterative Mode Coupling : With the heat source q 2 , a solution of the 

fundamental mode (1*2) is obtained which describes a simple Hadley cell and 
the rigid shell component of superrotation (Mayr and Harris, 1981 ). The 
meridonal velocity In that mode couples moment urn, through the Coriolis force, 
into the next higher mode (see the Equation (13) and the connection on the 
left hand side of the block diagram in Figure 6), and we obtain successive 
solutions for t>4,6... Choosing some limiting wave number, we then turn 
around in our iteration procedure and consider the feedback due to momentum 
coupling by the zonal velocities (see Equation (12) and the connections on the 
right hand side of the block diagram in Figure 6), until we return to the 
fundamental mode of superrotation. The results substantially confirm our 
earlier conclusions. This solution procedure shows in particular that the 
feedback from the zonal velocity component in the higher wave numbers (shorter 
wavelengths) contributes to accelerate the equatorial region in the lower 
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order inodes (longer wave lengths). Since the energy clearly comes from the 
heat source in the lowest order mode, our interpretation is that the feedback 
acts as a trapping mechanism. Under viscosity, velocity shears can 
numerically propagate in altitude, and the procedure eventually becomes 
unstable as one increases the number of iterations. 

(2) Simultaneous Mode Coupling; With a limited number of 12 modes, 
momentum coupling is considered in self consistent form by simultaneously 
solving a system of 60 second order differential equations (egs. 9 through 13 
for 1 - 2 to 24). The shortest wavelength covers about 15° in Latitude which 
is crude by comparison with the observed details in the Jovian circulation. 

It is again assumed that the lowest order heat source q 2 is driving the 
motions. 

Results are shown for the zonal velocity and the relative temperature 
amplitude in the form of latitudinal distributions at various altitudes 
(Figures 10 and 11). The meridional circulation, of the Ferrel-Thompson type, 
is Illustrated In Figure 12. One should expect that the dynamic properties of 
the atmosphere are sensitive to the tropospheric temperature lapse rates and 
the vertical and horizontal dif fusivltles as they vary with altitude. To 
understand the model fully and in relation to observations, will therefore 
require an extensive parametric analysis, the subject of a subsequent paper. 

In the lower troposphere the temperature lapse rate is almost 
adiabatic. The meridional circulation does not extend all the way down but 
instead Is broken up Into smaller cells (Figure 12) which is reflected in the 
complex structure of the zonal wind field (Figure 10). Higher up, where a 
larger temperature gradient is required, <5x10 ^ , to conduct the heat 
from the planetary Interior, the meridional circulation effectively fccuses 
and traps energy and angular momentum In the equatorial region. local 
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temperature maximum (Figure 10) and a well defined equatorial jet thus develop 
at the equator* Very high up in the atmosphere where the temperature lapse 
rate Is subadlabatlc. heat is diffusing downwards and effectively disperses 
the energy that Is transported upwards by the meridional circulation. As a 
result the latitudinal temperature distribution flatens and a blunt maximum 
develops at low latitudes where the solar energy Is deposited (Figure 11). 
Moreover, angular momentum Is carried to higher latitudes and broad maxima 
develop In the (eastward) zonal wind velocity. Due to upward convection In 
the form of thermal expansion from below and vertical as well as horizontal 
momentum diffusion, the equatorial region superrotates even at higher 
altitudes. 

In the Intermediate regions near 30 km, between the "dispersing" stable 
region above and the "focusing" unstable region below, the form of the zonal 
wind field resembles that of Jupiter. The horns near 10^ latitude are 
reproduced and secondary maxima develop near 30*^ and 60^. Lower down well 
within the unstable convection region at -20 km, the equatorial Jet Is larger 
by more than a factor of three and the form of the latitudinal velocity 
distribution resembles that of Saturn (Smith et al. , 1981). 

Based on this result we suggest that the large difference In magnitude 
between the equatorial jets of Jupiter and Saturn can In part be explained 
(Figure 13) by the temperature differences observed on both planets (Hanel et 
al., 1981). The clouds, which allow us to trace atmospheric motions, are 
formed on Jupiter near the tropopause where the transition occurs between two 
vastly different dynamic regimes one convectlvely stable the other one 
unstable. A complex Interfetence pattern thus develops In the zonal 
circulation. On Saturn which Is much colder, the clouds form at a much lower 
altitude, deep Inside the convection region (S^ < 0) where the meridional 
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circulation prevails In funnellng energy and moocntuffl Into the equatorial 
region. 

CONCLUSION 

The Important elements of the model are summarized in Figure 14. On 
Jupiter and Saturn global scale convection carries energy from the Interior 
•xkxkuueci wnere radiation to apace becomes Important. The 
tropospheric regions should therefore be convectlvely unstable with small 
superadlabatic Lapse rates on the order of 10~^deg/cm. Under this condition, 
the upward motions In the direct and Indirect atmospheric circulations Imposed 
by ln«oIatlon can then supply dynamic energy for "localized" equatorial 
heating which In turn drives equatorial superrotation. Preferential solar 
heating at low latitudes Induces the ordered structure In the atmospheric 
motions. The efficiency of this process Is determined by the degree of 
turbulence In the troposphere, the energy conversion. In the planetary Interior 
and the radiative loss from the upper troposphere and stratosphere. 

Differences in these global properties and the relative location of the clouds 
may contribute to the large differences between the zonal velocity fields on 
Jupiter and Saturn. 

Adopting a spectral model, the concept of mode coupling is used to 
describe differential rotation and equatorial superrotation In an atmosphere 
with zonal and meridional dscous stresses. Momentum Is thereby allowed to 
cascade from lower order modes, which "absorb" most of the solar energy, to 
higher order modes. At altitudes above the center of the meridional circula- 
tion, where the specific heat Input Is relatively large and the atmosphere Is 
convectlvely stable, angular momentum is effectively advected poleward, except 
near the equator where the meridional winds and the Coriolis force vanish and 
the horizontal momentum transport must be carried by horizontal diffusion 
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(viscous flow)* This causes the atmospheric rotation rate to decrease at low 
latitudes and increase at high latitudes relative to the rigid shell component 
of superrotation; inducing a tendency for the equatorial region to corotate* 
However* below the center of the meridional circulation and below the maximum 
energy input where the atmosphere tends to be convectively unstable* angular 
momentum (and energy) are advected from high to low latitudes, thus the 
equatorial region can be accelerated to the extent that equatorial Jets 
develop. Again* near the equator, angular momentum must be carried by 
horizontal diffusion; the equatorial minimum In the angular velocity on 
Jupiter (Figure 2) being evidence for down gradient diffusion toward the 
equator* This picture Is complicated by vertical energy and momentum 
diffusion as well as convection which can effectively extend the 
characteristic circulation signatures from the troposphere Into the 
stratosphere or visa versa. 

Our theoretical results suggest that four factors are Important In 
understanding the structure of the zonally symmetric circulation on Jupiter 
and Saturn: (1) Their large planetary size and rotation rate produce through 

mode/momentum-coupllng a multiple Ferre l-Thomson type circulation with upward 
motions at low latitudes. Differences In the gravitational acceleration as 
well as In the planetary radius and rotation rate of both planets 
significantly contribute to broaden the latitudinal structure of the zonal 
circulation on Saturn relative to Jupiter. (2) Energy from the Interior is 
diffusing radially outwards and maintains a superadlabatlc lapse rate* This 
provider the condition for trapping or focusing dynamic energy that Is 
advect.>d toward the equator In the lower leg of the meridional circulation. 

(3) The eddy dlffuslvlty Is sufficiently low to require large temperature and 
velocity gradients In the energy and momentum balances. A vertical eddy 
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diffusion coefficient between 10^ and 10^ may produce the conal yelocltles 
observed on Jupiter; but a amaller value is probably required to explain the 
zonal circulation on Saturn. (4) Since the clouds are the tracers of 
atmosphric motions, their location is very important. On Jupiter they occur 
near the tropopause, in the transition zone between convective ly stable and 
unstable regimes, and a complex interference pattern develops in the zonal 
circulation. In the atmosphere of Saturn which is colder, the clouds occur 
lower down well within the unstable convection region where the meridional 
circulation prevails in funneling energy and momentum into the equatorial 
region. This in turn is conducive Lo the development of n •q*»*torlal 

Jet. 

Mayr et al., 1981 argue that the large scale circulation on Jupiter 
provides the order in the formation of the small scale circulation features 
such as the great red spot and the white and brown ovals. Thus a relatively 
simple picture emerges for the Jovian circulation. The energy Jupiter 
radiates to space must be transported upwards through the troposphere. If 
that transport Is accelerated by the prevailing up%»rd motions in the solar 
driven meridlonlal circulation, eastward Jets develop such as observed in the 
equatorial region. But, if that vertical transport Is impeded by the pre- 
vailing downward motions in the solar driven multicellular meridional circula- 
tion, the atmosphere "reacts" and tends to maintain the process through the 
development of hurricanes. Dynamically induced by solar differential heating, 
a beautifully simple and ordered latitudinal structure, with alternating 
"stability" and "instability", is imposed on the troposphere, to form 
alternating zonal strata where hurrlcances In the form of red, bro«m, and 
white ovals are "forbidden" and "pennitted", respectively. 
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FIGURE CAPTIONS 


figurt 1 : Empirical model of the tonal wind velocitiae observed from Che 

Voyager spacecraft (C* Hur.c, private communication) in terms of vector 
spherical harmonics which are us«d as basis functions in the subsequent 
theoretical analysis. Data from both hemispheres are taken, but only the 
terms symmetrical with respect to the equator are presented (lotrer part) and 
considered for the syntheses. Syntheses up to wave number L ■ 12 and 34 are 
shown (upper port). The rigid shell component (t*2) , directly driven by 
solar differencial heatirg, dominates, and the lower order terms reproduce the 
equatorial Jet. Note the simple form of Che power spectrum; the Voyager I and 
II data being nearly identical. 

Figure 2 ; The zonal wind field based on vToyager and Earth based measurements 
is taken from the NASA publication Voyage to Jupiter (Morrison and Samz, 1980) 
and shows a minimum below /° latitude, indicative of downgradient diffusion 
toward the eqvutor. We consider an equatorial region (right side) in which 
the angular momentum budget is established with horizontal and vertical 
diffusion. Thus, a priori, there is no constraint on the direction or 
magnitude of Che zonal velocity at the equator and the question reduces to how 
angular momentum is effectively transported into the e:;uacorial region? 

Figure 3 ; Schematic Illustration of the mass (top)-and raomentuo (botcoQ)** 
budgets for the direct circulation producing rigid shell superrotation which 
is observed on *.*lrtually every planet in our «olar system. The balance is 
established between upward transport of momentum by convection (more being 
carried up%<ards at low latitudes than is carried do%mwards at high latitudes) 
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and diffusion. During "spin up" the atmosphere receives the angular momentum 
for superrotation from the planet or the planetary Inferior. 

Figure 4 ; Relative to the rigid shell component of superrotation (dashed 
11 angular momentum Is redistributed by the Indirect circulation vlthln 
the atmosphere. The momentum transport Is tomard the poles In the upper leg 
of themeiTldlonel cell thus Increasing the superrotation at mid latitudes (plus 
sign) and decreasing It near the equator (minus sign). In the lower leg of 
the merldlo!.'*' circulation this process Is reversed; angular momentum Is 
redistributed toward the equator and an equatorial jet can develop (plus .Ign) 
under certain conditions (See Figure 5). It Is thereby understood that In the 
equatorial region the horizontal momentum transport Is taken over by horizon- 
tal diffusion (Figure 2). 

Figure 5 : Ignoring the added complexity due to horizontal diffusion, the 

momentum balance between horizontal advectlon and vertical diffusion is 

schematically illustrated. The equator Is understood to be the equatorial 

region. Three different height distributions satisfy the budget. Case (a) Is 

Illustrated In Figure 4. In cases (b) and (c) the atmosphere superrotates and 

subrotates at all altlrudes. Tlie outcome depends on the energy budget 

Illustrated In the lower part. To sustain equatorial superrotation the 

meridional circulation must accumulate energy at the equator. The condition 

for this arises when the atmosohere Is convectively unstable, 

3To 

F) < 0, which Is presumably satisfied in the tropospheres of 
Jupiter and Saturn. 
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Figure 6 : Block diagram llluat racing Che mode /momenc urn coupling for Che 


apeccral model. The heaC aource In Che fundamenCal mode, q 2 t dlrecCly 
driving Che rigid shell componenc of auperroCaClon. Mode coupling Chrough Che 
solenoldal (meridional) veloclCy field In turn la driving Che Indirect circu- 
lation by cascading momencum from lower to higher order modes (left hand side 
of the chain link). A feedback Chrough the solenoldal or zonal velocity field 
(right side of Che chain link) tends Co Crap angular moroenCura in the lower 
order modes. 

Figure 7 ; Heuristic discussion of various scenarios for the zonal velocity 
field In Che equatorial region based on the equations (22) through (24). An 
equatorial heat source Is forcing upward motions and equatorward and poleward 
winds in Che lower and upper legs of the meridional circulation, respectively. 
In the direct circulation a rigid shell component of superrotatlon is produced 
(a). Superimposed on Chat fundamental mode the Indirect circulation produces 
differential rotation with zonal velocities at the equator Illustrated in the 
upper part of Figure 7. The sum of Che rigid shell component and differential 
rotation Is shown In the lower part of Figure 7. Cases (b) and (c) represent 
conditions where the atmosphere is convectlvely stable and unstable, 
respectively. Under case (d) the atmosphere Is unstable and stable In the 
lower and upper legs of the meridional circulation, respectively. 

Figure 8 ; The analytical model (30) through (35) Is used to describe the 
latitudinal distribution In the upper and lower legs of the circulation 
corresponding to case (d) In Figure 7. The powei spectra are shown in the 
lower part; for this display the amplitudes of the Indirect circulation are 
scaled by a factor of three. Considering (25) the solution represents local 
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conditions where the two altitude regimes ere wtsekly coupled. When coupling 
Is Important near the transition from one to the other regime a complex 
Interference pattern can develop such as observed on Jupiter where besides the 
equatorial Jet eastward Jets are also observed at higher latitudes. 

Figure 9 ; Input data for the average temperature, T^, and pressure, p^, the 
vertical eddy diffusion coefficient, k, and the heat Input, q 2 « In the 
troposphere near the 1 bar pressure level, the (negative) stability Is assumed 
to be Sq " **5xl0~^; Its absolute value decreasing toward lower altitudes so 
that the heat flux Is conserved. 

Figure 10 ; Latitudinal velocity distributions at various altitudes obtained 
from a numerical solution of the theoretical model. Note the large equatorial 
Jet at -20 km and the velocity maximum near 30^ at higher altitudes. In 
between a complex velocity field develops which has some resemblance with that 
of Jupiter. 

Figure 11 ; Same as Figure 10 except that the relative temperature amplitude 
Ij shown. Note the sharp maximum at the equator at altitudes where the narrow 
eastward Jet develops. At higher altitudes the maximum Is very broad. Indica- 
tive of poleward energy transport (adiabatic cooling) by the indirect circula- 
tion. 

Figure 12 ; The meridional circulation, based on out theoretical results, is 
Illustrated. We refer to it as a multicellular Ferre 1-Thonson circulation. 
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Figure 13 ; Based oa the results from Figure 10 » we suggest that the observed 
differences In temperature between Jupiter and Saturn (Hanel et al., 1981) are 
In part responsible for the large differences between the velocity fields 
observed on both planets (Smith et al., 1981)« On Jupiter, the clouds 
condense near the tropopause which Is In a transition region between the 
stable and unstable convection regimes above and below, where poleward and 
equatorward transport dominate, respectively. Thus a complex Interference 
patten develops such as seen In Figure 10 at 20 km. On Satun vdilch Is much 
colder, the clouds condense about SO km lower well within the unstable 
convection regime. The lotier leg of the meridional circulation thus dominates 
In funnellng energy and momentum Im > '..i. equatorial region, much like 
hurricanes are formed on a much smaller scale. This process is conducive Co 
Che development of a relatively large equatorial jet such as seen In Figure 10 
at -30 km which Is analogous to the conditions on Saturn. 

Figure 14 ; The important elements of the model are summarized In schematic 
form. On average the troposphere is convectlvely unstable; a superadlabatlc 
temperature lapse rate accomodates Che upward transport of heat from the 
planetary Interior which Is radiated to space at higher altitudes. Above Che 
tropopause Che solar heat Input dominates and energy Is conducted to lower 
altitudes where colllslonal excitation permits reradlatlon. Against this 
"homogeneous" background, the atmosphere is preferentially heated by solar 
radiation at low latitudes which brings order Into the large scale atmospheric 
motions, effectively polarizing them. Air Is rising at the equator and falls 
at higher latitudes. Thus angular momentum from the vast reservoir of the 
planet Is stored In the atmosphere during "spin up" causing It to superrotate 
on average (rigid shell component of superrotatlon). With the large rotation 
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rate and radius of tha planet tha Hadlay circulation is brokan up Into a 
nultlcallular maridional circulation with altarnatinp rising and falling 
motions in tha troposphere* Under conditions which are convectlvaly unstable) 
the upward motions, focused by the Coriolis force into narrow latitude regimes 
effectively funnel energy and momentum into the equatorial region, for 
example, which is conducive to the formation of an equatorial Jet* 
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